Transient Turbulence in Taylor-Couette Flow 
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Recent studies have brought into question the view that at sufficiently high Reynolds number tur- 
bulence is an asymptotic state. We present the first direct observation of the decay of turbulent states 
in Taylor-Couette flow with lifetimes spanning five orders of magnitude. We also show that there 
is a regime where Taylor-Couette flow shares many of the decay characteristics observed in other 
shear flows, including Poisson statistics and the coexistence of laminar and turbulent patches. Our 
data suggest that characteristic decay times increase super-exponentially with increasing Reynolds 
number but remain bounded in agreement with the most recent data from pipe flow and with a 
recent theoretical model. This suggests that, contrary to the prevailing view, turbulence in linearly 
stable shear flows may be generically transient. 

PACS numbers: 47.27.Cn, 47.20.Ft, 47.27.ed, 47.27.N- 



Given its wide applicability and importance, it is some- 
what surprising that there is still a rather incomplete un- 
derstanding of the physics underlying turbulence. While 
the statistical approach of the past century has provided 
great insight, it has proven to have limited quantitative 
predictive power, even for the simplest of flows. Dynam- 
ical systems has emerged as a complementary technique 
and has shown promise in clarifying the physics of tur- 
bulence In this picture, turbulent flows are viewed 
as trajectories in a high-dimensional phase space that 
wander between unstable solutions that coexist with the 
laminar solution [2|. 

Experiments in shear flows have shown that below the 
onset of turbulence, but at sufRciently high Reynolds 
number {Re = UL/v, where [/,L,andi/ are the veloc- 
ity and length scales of the problem and the kinematic 
viscosity of the fluid, respectively), states where well- 
defined regions of turbulent and laminar flow coexist can 
be reached by finite-amplitude perturbations to the lami- 
nar flow. For a range of Re these states decay back to the 
laminar state. The lifetime of an individual event is very 
sensitive to the details of the initial perturbation [1, |J. 
However, ensembles of experiments using similar pertur- 
bations have been found to have an exponential distri- 
bution of lifetimes for fixed Re [6| . This distribution 
is consistent with a chaotic repeller model 0, [l] ■ Some 
numerical simulations and experiments have concluded 
that at moderate Res the characteristic lifetimes of these 
states increase with increasing Re until at some critical 
Reynolds number Rec they diverge and turbulence be- 
comes sustained [1, 0, ■ Above Rec the dynamics of 
the turbulent state are typically associated with those of 
a chaotic attractor [lT|. Other experiments and simula- 
tions suggest that characteristic lifetimes do not diverge, 
instead increasing very rapidly but remaining finite at fi- 
nite Re [ll, [13, m, [l5| . This implies that turbulence is 
not a permanent state of the flow for any Re and is in- 



stead generically transient, if very long-lived. This result 
has important repercussions when one considers the pos- 
sibility of controlling the transition to/from turbulence 
by exploiting dynamical connections between the lami- 
nar and turbulent states [l6!| . The issue remains a subject 
of debate in part because the open flows (i.e., pipe flow 
and plane Couette flow) have been used to address this 
question. This limits the maximum time during which a 
turbulent episode may be observed and introduces com- 
plications associated with inlet conditions. Numerical 
exploration of this question is also problematic since it 
requires many highly-resolved simulations over long times 
to obtain significant statistics. 

We present the first measurements of lifetimes of tran- 
sitional states in Taylor-Couette flow (TCF). TCF (i.e., 
the flow between rotating, concentric cylinders) is most 
famous for a transition to turbulence preceded by a hier- 
archy of bifurcations that arise from centrifugal instabili- 




FIG. 1: Photographs of turbulent patches in TCF at Re = 
7500 with only the outer cylinder rotating. In this regime, 
turbulent patches coexist with the laminar flow and evolve in 
space and time. For all Re studied, these patches decay away 
in a probabilistic manner with a characteristic timescale r 
dependent on Re. The photographs show a 25 cm high region 
of the flow. 
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ties [13, Ell- However, TCF can bypass these instabilities 
and make a direct transition to turbulence. Several stud- 
ies have looked at transitional states that arise in certain 
regimes of counter-rotating TCF [l^ [2^ , but little work 
has been done when only the outer cylinder is allowed 
to rotate. We believe that this regime of TCF is ideally 
suited to the experimental investigation of the lifetimes 
of transitional states. 

In this regime, the transition to turbulence is abrupt 
and is characterized by spatially and temporally inter- 
mittent patches of turbulence that coexist with a laminar 
background (see Fig. [T]). Transition is observed to 
occur despite calculations that show that circular Cou- 
ette flow should be linearly stable for all outer cylinder 
rotation rates as long as the inner cylinder is stationary 
[2l| . These are features that are observed in other tran- 
sitional shear flows [13, [1^. However, because TCF is 
streamwise periodic, it allows for arbitrarily long obser- 
vation times, while avoiding the problems of contamina- 
tion from the inlet that plague plane Couette and pipe 
facilities. Furthermore, using this untested geometry al- 
lows us to test the generality of the trends observed by 
earlier studies of the decay times of transitional states. 

Our experimental apparatus was a vertical Taylor- 
Couette system with a glass outer cylinder of radius To 
— 7.620±0.008 cm and a brass inner cylinder of radius 
r,; — 6.638±0.001 cm. The gap width d — To — Vi was 
0.982 cm and the radius ratio 77 = ri/vo was 0.871. The 
cylinder ends were attached to the outer cylinder and the 
cross-stream aspect ratio T — L/d was 36. Each cylinder 
was driven by a stepper motor and their angular veloc- 
ities LUo^i could be set to better than 0.5%. The entire 
apparatus was submerged in a temperature bath set to 
20.00°C that varied by less than ±0.05°C over several 
days. 

Our working fluid was distilled water with 2.2% 
Kalliroscope AQ 1000 added by volume for flow visualiza- 
tion 24]. Adding Kalliroscope to water slightly changes 
its viscosity, so we measured it using a Cannon-Fenske 
routine viscometer. We found that it had a kinematic 
viscosity of 1.0298 ± 0.0025 mmVs at 20°C. We also 
studied the rheology of the suspension using a stress- 
controlled rheometer and found no deviations from New- 
tonian behavior. These results, along with the geometric 
characteristics of the apparatus, allowed us to set the 
outer cylinder Re to within 1%. Here, 

UL roUJo{ro~ri) 
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Like plane Couette and pipe flows, TCF (in the regime 
described above) requires a finite-amplitude perturbation 
to trigger transition. The laminar state was first pre- 
pared by accelerating the outer cylinder to the desired 
Re. Since at sufficiently high Re (~16,000, in our case) 
ambient noise is sufficient to cause spontaneous transi- 
tion, we limited our study to i?es significantly below this 
threshold. The system was allowed to run in Couette fiow 
for several minutes (~ 3 radial diffusion times td = d^ /v) 
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FIG. 2: The inner cylinder angular velocity as a function of 
time sliown above was measured using solid state rotation sen- 
sors. Experiments conducted at Re — 7647 showed measured 
lifetimes to be insensitive to variations in inner cylinder accel- 
eration, maximum velocity, and pulse duration so the profile 
above was used for experiments at all values of Re studied. 



in order to eliminate any transients. 

The flow was then perturbed by rapidly accelerating 
the inner cylinder in the direction opposite the rotation 
of the outer cylinder and immediately stopping it. Fig- 
ure [2] shows the angular velocity of the inner cylinder 
as a function of time. It is important to note that our 
perturbation is qualitatively different from those used in 
the pipe flow experiments [^ [T3 | and is similar to the 
quenching experiments conducted in plane Couette fiow 
y because it is not a small, localized perturbation but a 
large, global one where the whole flow is disturbed. After 
briefly exhibiting featureless turbulence, the flow relaxed 
to an intermittent state like the one shown in Fig. [1] 
The fiow was then observed until it relaminarized. 

A CCD camera captured video of the fiow at a reso- 
lution of 320 X 240 pixels. Except for the three highest 
Reynolds numbers studied, the video was captured at 30 
frames per second. At these Reynolds numbers the frame 
rate was reduced to 5 frames per second due to data stor- 
age limitations since individual events could last many 
hours. The low resolution of the video stream allowed it 
to be analyzed in real-time by a computer. 

First, the video stream was separated into ten frame 
segments and a frame from each segment was subtracted 
from a reference frame in the segment preceding it. The 
resulting image was thresholded to highlight only pix- 
els that were significantly different from the same pixel 
in the reference image to reduce noise. The elements of 
this binary image where added together and the resulting 
number N was compared to an empirically determined 
threshold, Nt- If N dropped below Nt for 30 seconds 
(turbulent patches were never observed to return after 
disappearing for more than a few of seconds), the system 
stopped acquiring data and prepared for the next run. 
The video data was saved and the accuracy of the auto- 
mated lifetime measurements was verified by hand for a 
random selection of the experimental trials. 
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FIG. 3: The probability of observing events whose lifetime is 
longer than time t decreases exponentially for long times with 
a decay constant that increases with increasing Re. Times are 
expressed in dimensionless time units tnd = d/roOJo- Only a 
fraction of the total events used to obtain the results in Fig. |4] 
are shown for clarity. The error bars indicate sampling error 
and are shown only at representative points. The solid lines 
show the weighted least squares linear fits to the tails of the 
distributions. 
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FIG. 4: Characteristic lifetimes as a function of Re in di- 
mensionless time units. * indicates data sets for which all 
the decays were directly observed. The points indicated with 
• correspond to data sets for which a maximum observation 
time was set as indicated in the text. Horizontal error bars 
represent absolute error in Re based on experimental limi- 
tations. Errors in logio t"^ are on the order of the symbol 
size and were estimated by the uncertainty calculated for the 
slope of the linear fits shown in Fig. [3] Inset: A plot of 
vs. Re on a linear scale shows that r^^ only as Re oo. 



After repeating the experiment many times (between 
200 and 1200 times depending on Re) at fixed Re, we cal- 
culated the cumulative probability distributions P(t, Re) 
shown in Fig. [31 P{t, Re) is the probability of observing 
turbulence at some time t after the initial perturbation 
for a fixed value of Re. 

These curves have two salient features. First, the 
distributions have exponentially decaying tails (i.e., 
P{t, Re) ^ exp{—{t — to)/T{Re)), where to is the time as- 
sociated with the initial formation of the turbulent state 
and T{Re) is a decay constant that depends on Re). This 
indicates that for long times the decay of turbulence is a 
Poisson process, which is a hallmark of a chaotic repeller 
0, Q . The same behavior has been reported in all previ- 
ous studies of the decay of turbulence in plane Couette 
and pipe fiows [1, i, [13, [H, [H, Q [H . r was calculated 
by fitting a straight line to the tail of each distribution. 
As a preliminary test of the role played by the end walls 
in the decay, we checked 150 long-lived events at various 
Res to see if turbulent patches were last seen at the top, 
middle, or bottom of the cylinder. We found that each 
location was equally likely. A more thorough study is un- 
derway to check the sensitivity of decay times to varying 
end conditions (e.g., variations in F). 

Second, the curves show that a significant fraction (~ 
30%) of the experiments relaminarized immediately af- 
ter the perturbation, almost independently of Re. Be- 
cause our perturbation consists of globally disorganizing 
the fluid, one can think of it as moving the system away 
from the laminar attractor in a random direction in phase 
space. If the perturbation puts the system in the basin of 



attraction of the turbulent state, the system will remain 
turbulent. If the perturbation leaves the system in the 
basin of attraction of the laminar state, the turbulence 
decays immediately. While the short-lived events where 
not used in determining the decay constants, the dura- 
tion of these events allowed us to estimate to. This was 
determined to be ~ 120 nondimensional time units. 

We observed that the characteristic lifetimes increased 
very rapidly with Re, but did not find a Re for which 
we could not observe decay. At Re = 9,922, we observed 
an event that lasted 29 hours before relaminarizing. The 
lifetimes became so long that we imposed maximum ob- 
servation times for the highest Reynolds numbers stud- 
ied. Because the experiments at lower Reynolds numbers 
were conducted first, the maximum observation times 
were chosen by linearly extrapolating from the decay con- 
stants calculated for the previous three experiments on 
a plot of logio r(i?e) vs. Re and multiplying this time 
by 1.5. This prevented us from directly observing the 
decays of the longest-lived events, but still left enough 
events of intermediate length to get good estimates of 
the characteristic lifetimes. To validate this procedure, 
we used it to measure the characteristic lifetime for Re = 
7,111. The resulting lifetime differed from that measured 
by observing all events by less than 1%. 

We also conducted five experimental runs of 300 events 
each at Re = 7647 to check the sensitivity of our results 
to the details of the perturbation. Increases in the pulse 
duration of up to a factor of 10 did not change our results 
significantly. Neither did decreasing the inner cylinder 
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acceleration by a factor of 60. The measured lifetimes 
were also insensitive to the maximum inner cylinder ve- 
locity as long as this was on the order of the speed of 
the outer cylinder. However, if the maximum speed was 
much lower than that, the fraction of events that relam- 
inarized immediately increase dramatically. This agrees 
with recent experiments in pipe flow |25j that showed 
that as long as the perturbation is large enough to cause 
the transition to turbulence the observed decay times do 
not depend on the details of the perturbation. 

Figure 2] shows the characteristic lifetimes r(i?e) as a 
function of Re. We fit the data with the various func- 
tional forms suggested by Hof et al. and found that 



r(i?e) ^ = exp{~exp{ci Re + C2)), 



(2) 



with ci = 3.61 X 10"'' and C2 = —0.59 best captured 
the trend (i.e., had the smallest residuals). While Eq. 
[2] seems to fit the observed trend, goodness-of-fit statis- 
tics indicate that it is not statistically significant fit (i.e., 
xt 3> 1). Therefore, we do not claim that it represents 
the actual functional dependence of r on Re, but only 
that lifetimes grow faster than exponentially but remain 
bounded. However, it should be noted that the scaling 
of Eq. [2] agrees with the most recent results for pipe flow 



14 1 and with the only theoretical prediction of T{Re) 
26|. As shown in Fig. Hj alternative fits to the data are 

possible (e.g., r^^ = exp{—{Re/ c^Y*^) with C5 = 3305 

and cg = 2.62) and differentiating between them requires 

many more decades of data. 

We have shown that when only the outer cylinder is 

allowed to rotate, transitional states in Taylor- Couette 



flow share many of the decay characteristics observed in 
other canonical shear flows. We showed that the decay of 
these states is a Poisson process, which suggests that the 
correct model for turbulence in this regime is a chaotic 
repeller. Furthermore, we have shown through direct ob- 
servation that characteristic decay times increase super- 
exponentially with increasing Re but remain bounded, in 
agreement with the most recent data from pipe flow 11411 
and with a recent theoretical prediciton for scaling [26l |. 
This suggests that contrary to the conventional view, tur- 
bulence may be generically transient for all Re in linearly 
stable shear flows. 

However, several open questions remain. As has been 
pointed out recently, it may be possible that turbulence 
can become sustained globally even if it is transient lo- 
cally in a manner similar to directed percolation [27j . 
This prediction comes from by a highly simplified model 
of plane Couette flow and must be verified experimen- 
tally for real fluid flows since it would be very computa- 
tionally expensive to address in direct numerical simula- 
tions. Other researchers have pointed out that the lami- 
nar state of flnite-sized Taylor-Couette system is nontriv- 
ial [12] and the role played by end-effects (e.g., Ekman 
pumping) in relaminarization is still unknown. Finally, 
it is unclear what effects curvature and rotation have on 
the decay of turbulence. Work is currently underway to 
address these issues. 
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